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Abstract 

We present the results on the source spectrum function for primary nuclei 
in galactic cosmic rays, where two distinct energy dependences are used for the 
source spectra. We discuss the evolution of the goodness of fit to B/C data with 
the propagation parameters and also show that the results are not much affected 
by a different choice for the diffusion scheme. We apply the constraints on the 
diffusion scheme as derived from stable nuclei to calculate the propagation of 
beta-radioactive isotopes. The diffusion model is refined to properly take into 
account the effect of a local bubble in the interstellar medium. Our calculations 
are compared to existing data, which prefer the local bubble description instead 
of the homogeneus scheme. 

1. Introduction 

Galactic cosmic rays detected with energies from 100 MeV/nuc to 100 
GeV/nuc were most probably produced by the acceleration of a low energy galactic 
population of nuclei and diffused by the turbulent magnetic field. The acceleration 
and diffusion processes have a magnetic origin, so that they should depend on 
rigidity. The rigidity function for the diffusion coefficient is given by quasi-linear 
theory as 



where the parameters K and 6 should ideally be given by the small-scale struc- 
ture of the magnetic field responsible for the diffusion. The spectrum just after 
acceleration and for a species j, Q^(TZ), depends on the details of the acceleration 
process, which are not clearly understood. However, a power-law seems to be 
preferred [1] and we assumed the following two possible forms: 
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dQ/dp = Fr a (L = 6 kpc) 
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Fig. 1. From top to bottom: for each best x 2 m the plane 5 — 7 (L = 6 kpc), the 
corresponding values of log(A'o), V c and V a are plotted. 

where the value of a is still debated (grossly from about 1.5 to 2.5). 
We propagate accelerated particles in a two-zones diffusion model, which has 
been described at length in [1, 2]. Briefly, our model takes into account spatial 
diffusion and convection wind {V c ), both acting in the thin matter disc and in the 
diffusive halo (this last has unknown size L kpc). Spallations over the interstellar 
matter, reacceleration (V^) and standard energy losses are set in the thin disc. 
The complete diffusion equation and its solutions may be found in [1] (see Eq. 4). 

2. Analysis of stable nuclei 

We tested our diffusion model by using data on the B/C ratio. The analysis 
was performed with six free parameters: Kq, 5,Vc,V a , L,a and with the imple- 
mentation of both the source spectra alternatively. We produced all the chain 
of primary and secondary nuclei by strating from Sulfur. Our aim is to compare 
predictions for the B/C spectrum with existing measurements, in particular with 
the 26 data points from heao3 [3]. 

In Fig. El we show the preferred values of the three diffusion parameters Kq, V c 
and V a , for each best x 2 m the 5 — 7 plane. Here we defined 7 = 5 + a. L has 
been fixed to 6 kpc but we checked that the behavior does not particulary depend 
on L. The source spectrum is the pure power law one of Eq. (j2J). The two upper 
panels show that the evolution of a does not affect Kq. On the other hand, we 
clearly see the (anti) correlation between K Q and 5 entering the diffusion coefficient 
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formula, because they need to give about the same normalization at high energy 
(K x .Ethresh ~ cte) . The lower left panel shows the values for the convective ve- 
locity Only very few configurations include V c = km s _1 , always when 5 = 0.3 
(for both types of source spectra). Increasing 7 and 5 at the same time makes 
V c change its trend. For small diffusion slope, convection is unfavored, as found 
in [4]. Finally, the Alfven velocity (lower right panel) V a doubles from 5 = 1.0 
to 0.3, whereas it is almost unchanged by a variation in the parameter 7. We 
found in [1] that the three parameters Kq, V c and V a behave very similarly with 
respect to a change in the source spectrum from "pure power law" to "modified" . 
In fact, its influence on the primary and secondary fluxes can be factored out if 
energy changes are discarded. The currently available data on B/C do not allow 
to discriminate clearly between these two shapes for the acceleration spectrum. 

3. Radioactive isotopes in the Local Bubble 

We studied the compatibility of our diffusion model with current data on 
/3-radioactive isotopes. These species diffuse on a typical distance / ra( j = \/Kjt 
before decaying. In this expression, not only the diffusion coefficient K, but 
also the lifetime 7r , depends on energy, due to the relativistic time stretch. In 
Table IT1 we give the values of some Z ra d for 10 Be, 26 Al and 36 C1. These species 





r (Myr) 


1 GeV/nuc 


10 GeV/nuc 


10 Be 


2.17 


220 pc 


950 pc 


26 Al 


1.31 


110 pc 


470 pc 


36 C1 


0.443 


56 pc 


250 pc 



Table 1. Rest frame lifetimes and corresponding values of l ra d for some (3 radioactive 
nuclei at two different energies ( Kq = 0.033 kpc 2 Myr" 1 and 5 = 0.6). 



are therefore very sensitive to the characteristics of the local interstellar medium 
(lism). The Solar System is embedded in an underdense region, usually called the 
Local Bubble (see [5] for references). The bubble leads substantially to a decrease 
in the spallation source term of the radioactive species. We modelled the bubble 
as a hole in the thin disc approximation and the radius of this hole is considered as 
an unknown parameter in the analysis. In [5] we provided the theoretical tools to 
treat the presence of the hole in the galactic disc in our two-zone diffusion model. 
The major result we found is that at the center of the bubble, the radioactive 
fluxes are decreased as 

jy/" r holc 

— — jr oc exp(-r hole // rad ) . 



Using the diffusion parameters allowed by the B/C data [2] to compute the 
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10 Be/ 9 Be and 36 C1/C1 ratios, we find that each of the radioactive nuclei points 
towards a bubble of radius < 100 pc, in relatively good agreement with direct 
observations. If these nuclei are considered simultaneously, only models with a 
bubble radius rhoie ~ 60 — 100 pc are consistent with the data and the case for 
r hoic = pc is disfavored. This is shown in Fig. |23 which is a projection of the 
parameter subspace allowed by B/C, 10 Be/ 9 Be and 36 C1/C1 on the L — 5 plane 
(left panel, no hole) or r ho i c — 5 plane (right panel, hole 7\ i e )- When the ratio 
26 A1/ 27 A1 is added to the analysis, the results become less clear [5], and it is sus- 
pected that the data (nuclear or astrophysical) on which they rely should not be 
trusted. 
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Fig. 2. Representation of the models compatible with B/C plus both 10 Be/ 9 Be and 
36 C1/C1 ACE 3-cj (open circles) and 1-a (filled circles) [7]. Left panel displays homo- 
geneous models (r^ie = 0) in the plane L — 5. Right panel displays inhomogeneous 
models (r^ie > 0) in the plane rhoic — S. 
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